Experimental Generation and Observation of Intrinsic Localized Spin Wave Modes in 

an Antiferromagnet 
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By driving with a microwave pulse the lowest frequency antiferromagnetic resonance of the quasi 1- 
D biaxial antiferromagnet (C2H5NH3)2CuCi4 into an unstable region intrinsic localized spin waves 
have been generated and detected in the spin wave gap. These findings are consistent with the 
prediction that nonlinearity plus lattice discreteness can lead to localized excitations with dimensions 
comparable to the lattice constant. 

PACS numbers: 05.45.-a, 76.50.+g,75.50.Ee 

Although solitons continue to play an important role in 
condensed matter physics in the last decade it was a ) 
recognized that nonlinearity plus lattice discreteness can 
lead to a different class of excitations with dimensions 
comparable to the lattice constant Such intrinsic 

localized modes (ILMs) have been identified in a vari- 
ety of classical molecular dynamics simulations m- 1] and 
with macroscopic mechanical [jl^] and electrical pi | mod- 
els, all of which ignore the possible role of quantum me- 
chanics 1 12 . Some effort has gone into identifying specific 
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condensed matter signatures as evidence of ILM produc- 
tion but all require intricate arguments: these include far 
infrared absorption [^3) , radiation ionization tracks |uj , 
the temperature dependent Mossbauer effect || and res- 
onant Raman scattering j^]. As yet evidence of exter- 
nally generated ILMs in a lattice of atomic dimensions is 
missing. The large amplitude modulational instability of 
an antiferromagnetic resonance (AFMR) for some anti- 
ferromagnets has been suggested fll|[l7j as a mechanism 
for the generation of intrinsic localized spin wave modes 
(ILSMs) and in this letter we describe the experimental 
observation and control of such nanoscale excitations. 

Because the driving field necessary to create ILMs 
via the modulational instability scales with the fre- 
quency of the antiferromagnetic resonance, wafmr, 
(C2H5NH3)2CuCi4 with the lowest frequency resonance 
at wafmr = 1.5 GHz was chosen for this first study. This 
antiferromagnet, often referred to as C(2)CuCLi, is a lay- 
ered organic material with a strong ferromagnetic cou- 
pling of the magnetic Cu 2+ ions within a layer and with 
a weak antiferromagnetic coupling between these layers. 
Because of this weak interlayer coupling the total spin in 
each layer can be represented by a classical one with re- 
spect to describing the lowest frequency mode dynamics. 
This spin system with its biaxial anisotropy is described 
in more detail in Ref. |l7[ . 

The C(2)CuCl4 single crystals were grown from aque- 
ous solution of ethylamine hydrochloride and copper (II) 
chloride in a closed vessel by slowly decreasing the tem- 
perature 18 1. For the measurement, platelets with well 
defined surfaces and a typical dimension of 3 x 3 x 0.5 mm 3 
were chosen. 



□ 



single loop, 
sample 



switch 




amplifier 




spectrum 
analyzer 



b) 



CD 

>, 
o 
c 

CT 
03 



1.4 







I I 1 1 

Ma 








! y F 








f \ 620 (is 


6 ms 






1 20 |is 

i i 





-2 







8 



2 4 6 

time (ms) 

FIG. 1. (a) Schematic experimental setup. Details are 
given in the text, (b) Experimental measurement road map. 
The horizontal line for t < identifies the AFMR frequency. 
During the high power driving pulse (shaded area) the detec- 
tion electronic is blocked. The solid lines A ... F show the 
AFMR for increasing pump-pulse power levels. The vertical 
dotted lines identify different time-cuts described in the text. 



The experimental setup is shown schematically in Fig. 
1(a). The first oscillator provides the high power pulse 
and the second oscillator the probe beam. Because the 
expected experimental signature of ILSM generation af- 
ter a strong microwave pulse near the wafmr is the 
breakup of the AFMR into a broad band, the genera- 
tion/detection system is tailored to produce a large AC 
field at the driving frequency cj e xcite and a high sensi- 
tivity over a broad band below but close to uj cxc itc- A 
single 3 mm diameter loop of copper wire is used as a 
non-resonant antenna for the excitation and pickup of 
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the broad band signal in reflection. A pump signal of up 
to 100 W can be obtained from the oscillator and solid 
state amplifier. To create short pulses a fast GaAs switch 
is used in front of the amplifier. The signal of a sec- 
ond, tunable oscillator is overlayed by means of a direc- 
tional coupler. With a second fast switch, synchronized 
by a digital delay generator, the reflected signal from the 
pump pulse is suppressed. The weak reflected signal from 
the second oscillator passes through two low noise ampli- 
fiers and is detected with a spectrum analyzer which was 
used as a variable bandwidth detector locked to the fre- 
quency of the probe oscillator. The minimum noise level 
of the system is — 120 dBm at 100 kHz bandwidth. In 
this measurement setup the absorption can be obtained 
as a function of time and frequency by first measuring the 
time dependent absorption at the frequency of the sec- 
ond oscillator, and then scanning this oscillator frequency 
during subsequent pump pulses. The time resolution is 
limited to 20 /is by the spectrum analyzer, and the sen- 
sitivity, by the flatness of the exciting loop impedance. 
Care was taken not to saturate any electronic circuits, as 
this could produce an artificial nonlinear response. The 
exciting coil and sample are immersed in pumped liquid 
helium at 1.2 K. 

To provide a road map to the experimental data-taking 
procedure a schematic view of the power and time de- 
pendence of the AFMR is shown in Fig. 1 (b). The 
AFMR frequency before the pulse, t < 0, is represented 
by the horizontal line. During the high power driving 
pulse (shaded area) the detection electronics is blocked. 
After the trailing edge (t — 0) of the driving pulse the 
spin wave is highly excited and its frequency is decreased 
due to the intrinsic nonlinearity of the spin Hamiltonian. 
From this state the system relaxes back to equilibrium 
with the longitudinal relaxation time Ti. The absorption 
by the spin system can be examined for different power 
levels of the driving pulse [lines labeled A ... F in Fig. 
1(b)] or for different time delays at a fixed power [dotted 
lines in Fig. 1(b)] or as a function of the driving pulse 
length. Measurements varying all of these experimen- 
tal parameters have been carried out and are described 
below. 

Figure 2 shows the absorption spectra 20 [is (3> T2 
of the AFMR) after the trailing end of a pump pulse 
for six different powers. The power of the 400 /is long 
driving pulse varied from 50 mW to 1.6 W. On the low 
frequency side of each spectrum are two magnetostatic 
volume modes while the shoulder on the high frequency 
side is a surface mode |lS|] . With increasing power lev- 
els the AFMR first collapses into a broad asymmetric 
shape and then at still higher powers returns to a sharp 
AFMR. Note that the magnetostatic modes do not show 
the same behavior as the AFMR. These experimental re- 
sults demonstrate that the AFMR is indeed unstable with 
increasing amplitude but only up to a specific amplitude 
while for still larger amplitudes the AFMR uniform mode 
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FIG. 2. Absorption spectra measured 20 /is after a 400 [is 
driving pulse for the different power levels. (A) 50 mW, (B) 
100 mW, (C) 160 mW, (D) 320 mW, (E) 500 mW, and (F) 
1.6 W of the 1.46 GHz (arrow) driving pulse. For clarity the 
spectra are shifted on the linear ordinate. The asymmetric 
wing on the low frequency side of the AFMR is the expected 
signature of localized mode production. The dot-dashed line 
marks the lowest frequency wafmr for an applied static field 
Hoc = 12 mT. 

again becomes stable. 

The results in Fig. 2 illustrate some of the criteria for 
the creation of ILSMs. First there is a minimum trans- 
verse amplitude above which the extended mode becomes 
unstable but then, surprisingly, there is also a maximum 
amplitude associated with this instability of the extended 
mode. This was not expected from our molecular dynam- 
ics simulations. Next the frequency interval of the insta- 
bility region is observed to become larger as the frequency 
shift Wexcite — wafmr becomes larger. For optimum con- 
ditions, localization could be observed for powers as low 
as 50 mW in a 400 /is driving pulse. Varying the pulse 
width from 50 [is to 400 [is while keeping the energy in 
the pulse fixed does not change these criteria. 

The narrow linewidth of curve F shown in Fig. 2 after 
the strongest driving pulse is a handy test to exclude 
any temperature related effect, such as the heating of the 
crystal by the intense microwave pulses. From low power 
linear measurements both the temperature dependence 
of the AFMR frequency, which goes to zero at the Neel 
temperature Tn = 10.2 K, and its linewidth are known. 
If the frequency shift of 75 MHz between the low power 
trace (dashed line in Fig. 2) and high power trace were 
caused by an increased temperature in the crystal, the 
line width should nearly double. This is not observed 
so temperature effects are not the source of the unusual 
lincshape results. 

The frequency response at several delay times for 
200 mW excitation in a 200 [is pulse is shown in Fig. 
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FIG. 3. Absorption spectra measured 20 /is, 220 [is, 420 /is, 
620 fis and 6 ms after a 200 fis long driving pulse at 1.46 GHz 
(arrow). The spectra are shifted on the ordinate by equal 
amounts. The dot-dashed line marks ojafmr as in Fig. 2. 



3. The five traces at 20 /is, 220 /is, 420 /is, 620 /is and 
6 ms show the evolution of this absorption feature from a 
broad resonance back to an AFMR at times smaller than 
Ti = 1.50 ms. For times larger than 2Ti only the AFMR 
survives. The radical difference between the line shapes 
at t = 20 /xs and t = 420 /is is consistent with the idea 
that the breakup into ILSMs is only possible when the 
effects of nonlinearity and dispersion are much stronger 
than the dissipation effect jl(| . 

MD simulations can be used to demonstrate that the 
extended spin wave is unstable against breakup into lo- 
calized modes for the specific amplitude of the extended 
wave produced in our experiments. Details on molecular 
dynamics procedures for simulating localized spin waves 
in C(2)CuCl 4 can be found in Ref. @. The solid line 
in Fig. 4(a) shows the calculated frequency dependence 
of the AFMR as a function of its transverse amplitude 
S y in the hard axis direction. To compare this simu- 
lation to our experiment, we extract from the shift of 
75 MHz or 0.05wafmr of trace (F) in Fig. 2 the spin 
wave amplitude S y = 0.06 as marked with (F) in Fig. 
4(a). Correspondingly the points A to F superimposed 
on this curve identify different AFMR frequency shifts 
observed in the high power measurements of Fig. 2. To 
test for the instability threshold the time dependent evo- 
lution of the energy density for a simulation of a 250 
spin antiferromagnetic chain is calculated, starting from 
an extended wave with fixed transverse amplitude plus 
random noise (5S n ) — 0.0025. For the whole range (A) to 
(F) covered by the experiment the extended spin wave is 
unstable and breaks up into localized modes with widths 
extending over about 10 lattice constants. Thus for these 
microwave powers the AFMR is driven into the interest- 
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FIG. 4. (a) Calculated frequency of the uniform spin wave 
mode as function of its transverse amplitude. (A) to (F) mark 
the experimentally determined position of the extended mode 
after the driving pulse from the corresponding traces in Fig. 
2. (b) AFMR plus ILSM absorption spectra determined from 
molecular dynamics simulations corresponding to three dif- 
ferent microwave power levels. Starting from extended spin 
waves with a given transverse amplitude of Sy = 0.01, 0.03 
and 0.06 (arrow positions a, j3, 7 in Fig. 4a) the uniform mode 
breaks up into ILSMs within a time of 200tafmr- From the 
evolution over the subsequent time interval At — 800tafmr 
the spectra are obtained via the auto-correlation of the net 
transverse magnetic moment as described in the text. 



ing nonlinear region. 

To compare the absorption spectrum A(oj) cx wx" 
measured in the experiment with results from molecular 
dynamic simulations, the imaginary part of the dynamic 
magnetic susceptibility is calculated using the Kubo ex- 
pression pc| where the absorption is proportional to the 
Fourier transform of the auto-correlation function of the 
net magnetic moment, M y (t): 

/•OO 

X y{uj) oc uj / dt{M y (t + t)M y {t)) exp (iut ) . (1) 



For our special case the simulations are started with 
an extended wave of a given amplitude plus some small 
amount of random noise. The extended wave is unstable 
against breakup into localized excitations and transforms 
after approximately 100 periods of the antiferromagnetic 
resonance into a broad spectrum of ILMs. From 200 to 
IOOOtafmr the evolution of the net magnetic moment is 
recorded and then via Eq. 1 the imaginary part of the 
dynamic magnetic susceptibility is found. 

The resulting calculated absorption spectra for a chain 
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of 1000 spins at three different powers are shown in Fig. 
4(b). To generate these spectra the starting transverse 
amplitudes of the AFMR are Sy = 0.01, 0.03, and 0.06, 
respectively. Each curve is averaged over 6 simulation 
runs to remove arbitrary spikes associated with single 
long-living ILSMs. The broad wing of the asymmetric 
band is due to the statistical distribution of the ILMs 
with different degrees of localization and thus different 
frequencies. To compare the simulations on the long 
chain with the spectra obtained in the experiment, the 
frequency axis in Fig. 4(b) was scaled to the same width 
relative to wafmr as given in Figs. 2 and 3. Both the 
shape and width of the experimental onset results are in 
good agreement with theoretical simulations. Moreover 
the absorption spectra shown in Fig. 4(b) are obtained 
starting from amplitudes of the extended wave which are 
in the right range for traces A to F. 

The disappearance of the instability at large ampli- 
tudes is an experimental feature not represented by curve 
7 in these model calculations. Since the values of the 
relaxation times taken from the experiment are Ti = 
1.50 ms and T2 « 0.1 /is, with the latter depending on the 
surface quality of the samples, T2 is much smaller than 
the pulse widths used in our experiments so that energy 
is transferred to other degenerate spin waves during the 
microwave pulse. Our experimental results indicate that 
at the powers where the uniform mode instability occurs 
this transfer docs not influence the positive curvature of 
the dispersion curve which is required for the instabilty. 
At the highest powers shown in Fig. 2 the uniform mode 
again becomes stable indicating that the finite wavevec- 
tor spin waves have become so heavily populated during 
the microwave pulse that the dispersion curve now has 
negative curvature at the uniform mode, a condition for 
stability of the mode at large amplitude (h| . Because it 
takes a finite time for the uniform mode to break up into 
ILSMs both features can appear at intermediate pow- 
ers such as displayed in trace D in Fig. 2. Here ILSMs 
are formed during the early part of the pulse while the 
population in the degenerate modes is still small, and 
they remain isolated in frequency space above the disper- 
sion curve after the degenerate mode population becomes 
large so both features are seen. 

In this series of microwave experiments the instability 
that appears, when the lowest lying AFMR of C(2)CuCl4 
is driven to larger amplitudes, has been used to gener- 
ate nonlinear excitations which are localized on a nano- 
lcngth scale. To distinguish between the ILSMs of clas- 
sical simulations and the excitations observed in experi- 
ment the latter will be identified as 'anons'. The hallmark 
experimental feature of the uniform mode breakup into 
anons is a broad and asymmetric spectral band below 
the AFMR. Classical MD simulations on a long 1-D an- 
tiferromagnetic spin chain show that both the amplitude 
of the extended spin wave at which it becomes unstable 
to breakup into localized modes and the resulting spec- 



tral shape of the ILSMs are in good agreement with the 
microwave results. At still higher experimental powers 
the AFMR is observed to become stable again, a feature 
related to the fact that the pump pulse is longer than T2. 

Simulations have played an important role in identi- 
fying the most economical experimental pathway for the 
detection of anons in real solids. This interplay may be 
expected to continue since additional MD studies on a 
C(2)CuCl4 spin chain in the presence of a magnetic field 
gradient indicate directed anon transport should be pos- 
sible. 
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